Singularities on the brane aren't isotropic 



Marco Bruni 1 and Peter K. S. Dunsby 2 ' 3 
1 Institute of Cosmology and Gravitation, University of Portsmouth, Mercantile House, Portsmouth POl 2EG, U.K. 
2 Department of Mathematics and Applied Mathematics, University of Cape Town, Rondebosh 7701, Cape Town, South Africa. 
3 South African Astronomical Observatory, Observatory 7925, Cape Town, South Africa. 

(July 20, 2002) 

Recent studies of homogeneous anisotropic universe models in the brane world scenario show that 
the cosmological singularity in this context is isotropic. It has therefore been suggested that this 
may be a generic feature of singularities on the brane, even in the inhomogeneous case. Using 
a perturbative approach, we show that this is not the case. As in the GR case, the presence of 
decaying modes in the perturbations signal the instability (in the past) of the isotropic singularity. 
The brane universe is therefore not born with isotropy built in: as in standard cosmology, the 
observed large-scale isotropy and homogeneity remains to be explained. 



I. INTRODUCTION 

The brane world scenario has recently received atten- 
tion as a possible string inspired cosmology (see |Q] for 
a review). In this scenario the observable universe is 
a 4-dimensional (4-d) slice, the brane, in a higher di- 
mensional spacetime, the bulk. Here we consider the 
particular implementation developed in in order to 
generalise a previous model by Randall and Sundrum 
P|, where the bulk is 5-dimensional (5-d) and contains 
only a cosmological constant, assumed to be negative. In 
this context various authors ^-^] have considered an ho- 
mogeneous and anisotropic brane, finding an intriguing 
result: unlike general relativity, where in general the cos- 
mological singularity is anisotropic, the past attractor for 
homogeneous anisotropic models in the brane is a sim- 
ple Robertson- Walker model. In particular, it is found 
in that this is also true for Bianchi IX models, as 
well as for some inhomogeneous models. 

In general relativity the Bclinski-Lifshitz-Kalatnikov 
conjecture [|l0| suggests that the Bianchi IX behaviour 
in the vicinity of the singularity is general, i.e. that the 
approach to the cosmological singularity in a generic in- 
homogeneous universe model should locally be the same 
as in Bianchi IX. 

Building on this and on the fact that Bianchi IX mod- 
els in the brane have isotropic singularities, it has been 
suggested in |||] that the isotropic singularity could be 
a generic feature of brane cosmological models. 

A well known problem of cosmology is to explain the 
very high degree of isotropy observed in the CMB. In a 
theory such as general relativity, where isotropy is a spe- 
cial rather than generic feature of cosmological models, 
we need a dynamical mechanism able to produce isotropy. 
Inflation was proposed, among other reasons, as a way 
to isotropise the universe. Inflation is effective in this 
sense, ^ but it needs homogeneous enough initial data 



* There are perturbative proofs of the cosmic no-hair conjec- 



in order to start at all flj. Although one can adopt 
the view that one such homogeneous enough patch in an 
otherwise inhomogenous initial universe is all we need to 
explain what we observe, this seems somehow unsatis- 
factory: the isotropy problem remains open in standard 
cosmology. 

If the conjecture in could be proved correct, brane 
cosmology would have the very attractive feature of hav- 
ing isotropy built in. Inflation in this context would still 
be the most likely way of producing the fluctuations seen 
in the CMB, but there would be no need of special initial 
conditions for it to startjj] Also, Penrose conjecture JlJ] 
on gravitational entropy and an initially vanishing Weyl 
tensor would be satisfied, cf. jllll . 

Unfortunately for the brane scenario, we prove here 
that the past attractor of homogeneous models fund 
in H|§ is unstable in the past to generic (i.e. anisotropic 
and inhomogeneous) perturbations. As in general relativ- 
ity, there is a decaying mode in scalar perturbations that 
grows unbounded in the past and that signal, in the con- 
text of linear perturbation theory, that anisotropy also 
grows unbounded as t — » 0. 

In order to prove this, we specifically look at pertur- 
bation modes of the dimensionless shear a/H, and show 
that there is a decaying mode (growing in the past) in this 
quantity. Ours is a large-scale analysis, at a time when 
physical scales of perturbations are much larger than the 
Hubble radius, A 3> H _1 (equivalent to neglecting Lapla- 
cian terms in the evolution equations). This may seem re- 
strictive, but this is not the case for the non-inflationary 
perfect fluid models tha are relevant to our discussion. 
Indeed in this case any wavelength A, smaller than H^ 1 



ture, i.e. classical perturbations in inflationary models with 
a scalar field or cosmological constant are swept out, as well 
as (local) proofs for homogeneous and inhomogeneous models 
(see e.g . pi| and references therein). 

^See |13| for CMB limits on the anisotropy in brane 
cosmology. 
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at a given time, becomes much larger than H~ x at ear- 
lier enough times. Because of this crucial property of 
perturbations for non-inflationary models our analysis is 
completely general, i.e. valid for any A as t — > 0. 

For the most part we follow the notation and conven- 
tion of [0,[l6|,[l7j . In section 2 we briefly summarise those 
results on general brane dynamics that are relevant for 
the following discussion (see Jl],||,|l6| for more details and 
other references). In section 3 we present the equations 
for scalar perturbations and derive the large-scale evolu- 
tion for the gauge-invariant density perturbation variable 
in the high energy limit. In particular we highlight the 
decaying mode that grows unbounded in the past. In 
section 4 we show that a corresponding mode in the di- 
mensionless shear a/H also grows in the past. In section 
5 we draw our conclusions. 



II. BRANE DYNAMICS 

The implementation of the branc-world scenario con- 
sidered in Q assumes that the whole spacetime is 5-d 
and governed by 5-d Einstein held equations: 



Gal 



-Ag AB + Sab(x) (-^9ab +T AB ) ■ (1) 



These represent a 4-d brane at \ — living in a bulk 
with metric c/ab and cosmological constant A; k 2 is the 
5-d gravitational constant, A is the brane tension, gAB 
and Tab are respectively the metric and the energy- 
momentum on the brane. The 4-d field equations induced 
on the brane are derived geometrically from (Q) assum- 
ing a Z2 symmetry with the brane as the fixed point, 
leading to modified Einstein equations with new terms 
representing bulk effects: 

G ' p, v = —Ag^ + K 2 T /Jt/ + k^S^v — £^ u , (2) 

where as usual k 2 = 8ir/M 2 . The various physical con- 
stants and parameters appearing in the equations above 
are not independent, but related to each other by 



A = 6, 



A 



\k 2 (A+ \k 2 \ 2 



(3) 



The tensor represents non-linear matter corrections 
given by 



Q 1 rp Oirp 1 rp rp(X 



3T aP T a P - (T a a ) 



(4) 



£^ v is the projection on the brane of the 5-d Weyl tensor; 
although the whole dynamics is 5-d and given by (Q), 
from the 4-d point of view is a non-local source term 
term that carries bulk effects onto the brane. If is the 
4- velocity of matter and h^ v = g^ + u^u v projects into 
the comoving rest space, one can decompose £^ u as KLEol: 



K A 



(5) 



effectively as if it was a trace-less energy momentum ten- 
sor with energy density U, energy flux Q M and anisotropic 
pressure V^v The brane energy- momentum tensor sep- 
arately satisfies the conservation equations, V v T flv = 0. 
Assuming a perfect fluid (or minimally coupled scalar 
field) we have the usual results: 



p + e{p+ P ) = 0, 

D liP +(p + p)A li = 0, 



(6) 
(7) 



where a dot denotes u v V v , 6 



D M u M is the volume ex- 



pansion, = u M is the 4-acceleration, and D M denotes 
the spatially projected covariant derivative. The con- 
tracted Bianchi identities on the brane then imply that 
the projected Weyl tensor £ M „ 
straints 



and Sftv obey the con- 



fix 2 

V O^u — — V D^v 



(8) 



which show how the non-local bulk effects are sourced 
by the evolution and spatial inhomogeneity of the brane 
matter content. Finally, using Eqs. (^)-(Q), these 
can be turned in propagation equations for the non- 
local energy density U and energy flux Q M . Neglecting 
terms quadratic in covariant variables that are gauge- 
invariant perturbations with respect to a Robertson- 
Walker isotropic background M,n8Ul9], these are:n 



\QU + D M <2„ = , 



D^ = -i K 4 (p + p)D M(0 , 



(9) 



(10) 



where H = a/a (= |0) is the Hubble expansion of the 
background. Using a standard decompositions of per- 
turbations, the anisotropic term V^ v carries in general 
contributions from scalar, vector and tensor modes; the 
latter however, satisfying YfV^v = 0, does not contribute 
the non-local conservation equations above. 

In the background, the Raychaudhuri equation is 



H 



-H 2 - — 
6 



p + 3p 

Ia- 



■3p) 



K 2 X 



-U n i± 



(11) 



* Strictly speaking, the variables denned in II8 M20I a nd those 
defined in the same way in the brane context ]l|, |llfpi| are 4-d; 
they are however easily generalised to 5-d. Bardeen like vari- 
ables have been defined in 5-d in order to carry out 
a brane-bulk analysis, e.g. see |2^], but the relation between 
these and the covariant ones used here has not yet been es- 
tablished; cf. |l9|j20| for this relation in 4-d general relativity. 
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where the solution for U follows from Eq. (^) (considering 
the zero order background only), a a is the initial scale 
factor and U = U{a ). The first integral of this equation 
is the generalized Friedmann equation on the brane (K = 
0,±1): 



K 2 X 



u [ a - 



(12) 



The high energy regime is defined by p 3> A. In this limit 
one obtains flat models dominated by the non-linear p 
term [EI, 



OA 



P 



(13) 



where we fix an arbitrary initial condition by choosing 
a Q = a{t Q ) = 1, and as usual w = p/p. This models are 
represented by a stationary (equilibrum) point, denoted 
J-f,, in the phase space of homogeneous Bianchi models 
consider in HH , as well as in the phase space of a special 
class of inhomogeneous Gi cosmological models. In both 
cases Tb is found to be the source, or past attractor, for 
the generic dynamics for w > (w — is also included 
in the homogeneous case) , consistently with . 

Finally, we note that the condition for inflation in gen- 
eral is fl29l 



1 2p + A 

w < -o — — r 
3 V P + A 



(14) 



As p/X — > oo this becomes w < — |, while the general 
relativity condition w < — i is recovered as p/X — > 0. 



III. DENSITY PERTURBATIONS 

Scalar gauge-invariant perturbations can be described 
using covariantly defined variables (see |l^,[l{| and ref- 
erences therein). In the brane scenario this formalism 
has been developed in Jl],|l6| (see also [pjj); here we shall 
follow the same approach, with minor modifications, and 
we refer to these papers for definitions. 

We can completely characterize scalar perturbations 
on the brane with four variables, A, C, U and Q, rep- 
resenting respectively the matter density perturbation, a 
convenient 3-curvature perturbation, the perturbation of 
the Weyl energy density, and a Weyl energy flux variable 
related to in (||) (all 3-quantities are defined with re- 
spect to a single 4- velocity field, that of matter The 
dynamics of these quantities is given by 



K 2 p(l + w) 



1 



+ 



C 



(1 + w) 
Aa 2 H 
Aa 2 Hc 2 
1 + w 



C 



2H 
p(l + w) 



2H 



P 
X 

6 
X^ 2 



U, 



(15) 



D 2 A 



12a 3 p\ , 



f 72a 2 Hc 2 \ / U 
V l + w J \Xk 2 



U= (3w-l)HU 
4c? 



l + w)\p) \3pJ 

(1 - 3w)HQ - —U - |aD 2 P 
3a 

1 lY 8c 2 \ U 4 , 
6a [ \ 1 + w J p 



A. 



(16) 



(17) 



(18) 



Note that this system of equation is homogeneous in the 
four chosen variables, except for the P term in (p^); 
this represent the contribution from the anisotropic Weyl 
stress, and since there is no evolution equation for P, in 
general one should either determine P from the full 5-d 
dynamics, or make some ansatz; otherwise, in general 
the system above is not closed and one cannot find so- 
lutions. Finally, the variable Z in ((T^), used in |fl^6 17 
and characterising the perturbation of the expansion, is 
related by A, U and C by the constraint 

C = 2n 2 a 2 p ( 1 + A + ^a 2 pU - 4a 2 HZ , (19) 
V XJ Xk z 

arising from the Gauss-Codazzi constraint in the brane. 

In the following we want to study the stability proper- 
ties of the models Tb, Eq. ([l3]), against generic inhomoge- 
neous and anisotropic perturbations, for values of w > 0. 
We see from (|l^) that these models are non-inflationary 
for w > 0. Thus, as we said in the introduction, for 
any A we only need to study the large-scale evolution of 
the variables above. We can either neglect the Lapla- 
cian terms in Eqs. (15)-(p"8|), or use a harmonic expan- 
sion (Fourier in our flat space case) and neglect terms 
H~ 2 /X 2 <C 1: the resulting equations are the same. 

Fortunately, in the large-scale limit the D 2 Q term in 
( |l7| ) is negligible, and one obtains a closed system for 
the density perturbations A and U and the curvature 
variable C . Besides, to the extent that the P contribution 
to (U|) is also negligible in this limit, cf. p^j, Q can also 
be determined p 16,171]. In addition, we now restric our 
analysis to the case U = 0: we will comment on the 
reliability of this assumption in the conclusions. 

Introducing the new variables U = U/(k 4 p) and Q = 
aHQ/ K A p and denoting with a prime the derivative with 
respect to dr = dln(a), we have 



A' = 



3w — 



K 2 P (1 



2H 2 



A 



(l + w) c 3(l + w)K 2 p 2 ~ r 



c 

U' 



4a 2 H 2 
0, 

2{3w 



XH 2 



1)U. 



(2-3w)Q-^U -^(l+w)A 



(20) 

(21) 
(22) 

(23) 



We see therefore that the evolution in the long wave- 
length limit, with U = 0, considerably simplifies: it is de- 
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scribed by a single first order equation for the density per- 
turbation and two first integrals arising from (|2l|)-(22). 
The first, C , represent the large-scale constant spatial 
curvature perturbation; the second, U = Ua~ 2 ( 1+3w \ 
is a first integral for the Weyl energy density pertur- 
bation. These source the density perturbation; there- 
fore, A has three modes, as opposed to the two arising 
in GR in the case of a single fluid. The first two are 
analogous to the GR ones: the first is the mode aris- 
ing from the C = U = initial condition, typically a 
"decaying mode", the second is the curvature adiabatic 
mode generated by Cq ; the third mode is a peculiarity of 
the brane scenario, is generated by Uq, and represent an 
isocurvature (or entropy) perturbation. Although this 
is not particularly clear from the treatment used here 
and in [p], p^|Jr^ ] , this mode is indeed due to the different 
intrinsic 4-velocities that the "Weyl fluid" and matter 
have in the perturbed spacetime (which is the cause of 
the presence of the energy flux variable Q). In complete 
analogy with the case of two fluids in GR, this then gcn- 
erates a "relative entropy perturbation" (see e.g. [ p0|j23| 
and [[l7j). Finally, and again as in GR, we remark that 
the adiabatic and isocurvature perturbations evolve in- 
dependently as described above only in the long wave- 
length limit: it is indeed clear from Eqs. (15)-([L8[) than 
in general the two modes are coupled and non-vanishing 
even starting from C Q = U a = 0. Conversely the de- 
caying mode should be seen as arising from special ini- 
tial conditions that lead to vanishing values of C and U 
when A 3> H . Usually this mode is decaying forward 
in time, and therefore neglected in structure formation 
studies, while the growing mode is the interesting one. 
In studying the question of homogeneity and isotropy of 
the brane at early times we are interested in running the 
equations above backwords in time, and it is the decay- 
ing mode that plays a crucial role for most values of w, 
as we are now going to show. 

Let us now restrict our analysis to the high-energy 
regime that dominates at very early times, when the 
background is given by the models Th, Eq. fll3|). Us- 
ing C and U , in the limit p > A the evolution of A is 
determined by the simple equation 

A' = -3A + j(l + wfC Q a q - 18(1 + w)U a r , (24) 

from which the three density perturbation modes are im- 
mediately determined. They are 

a a P 9(I + w)3 18(1 + w) r - r 

A = A a^ + - — — -C Q a q —=-U a a , (25) 

4 bw + 7 bw + 5 



where 



q = 6w + 4 



r = 2(1 + 3w) 



(26) 



g < ■£> iu < -| , r<Q^w<--, (27) 

and A is the constant of integration associated with the 
decaying mode. This shows that independently of the 



value of w, there is always a large-scale mode that grows 
unbounded in the past. 

The solution for Q can also be easily determined: 



6 5 — 3w 
3 (1 + w) 4 



M + ±( 2 + *?\w (28) 



8 (7 + 6w)(2 + 9w) 27w(5 + 6w) 

for and 

Q = Qoa 2 + ^A a- 3 - ^C Q a 4 + ^U Q a 2 \na (29) 

for w — 0. where 



s = 2 - 3w , s < ^ w > - , 



(30) 



and Q is a constant of integration representing the ho- 
mogeneous solution to equation (p3|). 



IV. THE EXPANSION NORMALISED SHEAR 

Like in the standard GR case, all the gauge-invariant 
geometric and kinematic quantities can be expressed in 
terms of A 28 in the context of linear perturbation the- 
ory. The key covariant variable related to the issue of 
isotropization in the past is the expansion normalised 
shear 0, 



S a 6 — 



Cafe 

IT 



(31) 



The scalar contribution to this quantity is obtained by 
taking its total spatial divergence pl| : 



a 2 D a D b a ab 
H 



(32) 



Using the shear constraint equation (equation (97) in 0] ) 
it is easy to show that S can be expressed in terms of A, 
C U and Q : 

S = 2A - -4t^ c + l2 U + 36( 3 , ( 33 ) 

and using (g5|) and ( gg ) together with the solutions for 
U and Q we obtain 



S = 



10 



5 — 3w 4 



3 (1 + w) 2 {99w 2 + 153w + 40) 



(6w + 7)(2 + 9w) 



8 18w 2 + 15w + 2 ~ _ ~ s . . 
3 w{b2 + 5) 



for w and 

E = ^A oa - 3 
5 



51 
28 



24 ~ 

C D a 4 + —U (1 + 2 In a) a 2 + 36Q G a 2 



(35) 

for w = 0. The presence of the decaying mode p = —3 
in ( |3"4| ) and ([55]) proves that E grows in the past. This 
completes our proof that the past attractor T b of homo- 
geneous models is unstable in the past against anisotropic 
and inhomogeneous perturbations. 
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V. CONCLUSIONS 

From a dynamical system point of view the past attrac- 
tor Th for brane homogeneous cosmological models found 
in [|]U is a fixed point in the phase space of these mod- 
els. This phase space may be though of as an invariant 
submanifold within an higher dimensional phase space 
for more general inhomogeneous models. The conjecture 
in d|| is equivalent to saying, in this dynamical system 
language, that Tb is the (local) past attractor for generic 
trajectories in this higher dimensional phase space. Our 
analyisis can be seen as an exploration of the neighbor- 
hood of Tb out of the invariant submanifold explored 
in We have found that Tb is unstable in the past 

to generic anisotropic and inhomogeneous perturbations 
of non-inflationary perfect fluid models with p = wp, for 
any value of w, using a large-scale A 3> H^ 1 approxi- 
mation that we have motivated in the introduction and 
is not restricting the validity of our analysis. The insta- 
bility of Tb we have found is fundamentally due, among 
other modes that may be stable or not depending on the 
value of w, to a decaying mode in the density pertur- 
bation that blows up in the past, A ~ a -3 , in a way 
independent of w and that, like in general relativity, is 
the signal of the unbounded growth of the dimensionless 
shear a/H as t — > 0, as proved in section 4. 

We have considered here only the case of vanishing 
background Weyl energy density, U = 0. This assump- 
tion considerably simplifies the analysis, but it is an easy 
guess that our results will remain true for W^O. Indeed 
when U J~b still remains the past attractor of the 
isotropic models whose stability in the past we want to 
examine, as is clear from the Friedmann equation (12). 
In other words, our analysis is restricted to the invari- 
ant submanifold U — of the larger phase space with 
U 7^ 0, but this submanifold is asymptotically stable 
against U perturbations. 

Finally, it has recently been suggested j3(J that the 
quantity Uq in ( |l2"| ) is only asymptotically a constant, 



while Uq 



at high enough energies. As it clear from 



( |l2| ) , even more in this case Tb still remains the relevant 
past attractor of isotropic models. We believe therefore 
that our analyisis should remain valid also in this case. A 
more complete analysis including this issue and therefore 
that of U ^ will be the subject of a future investigation. 
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